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Detection of B-mode Polarization in the Cosmic Microwave Background 
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Gravitational lensing of the cosmic microwave background generates a curl pattern in the observed 
polarization. This “B-mode” signal provides a measure of the projected mass distribution over the 
entire observable Universe and also acts as a contaminant for the measurement of primordial gravity- 
wave signals. In this Letter we present the first detection of gravitational lensing B modes, using 
first-season data from the polarization-sensitive receiver on the South Pole Telescope (SPT pol). We 
construct a template for the lensing B-mode signal by combining E-mode polarization measured by 
SPTpol with estimates of the lensing potential from a Herschel-SPIRE map of the cosmic infrared 
background. We compare this template to the B modes measured directly by SPTpol, finding a 
non-zero correlation at 7.70 significance. The correlation has an amplitude and scale-dependence 
consistent with theoretical expectations, is robust with respect to analysis choices, and constitutes 
the first measurement of a powerful cosmological observable. 
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Introduction: Maps of the cosmic microwave back- 
ground (CMB) [1] polarization anisotropies are naturally 
decomposed into curl-free E modes and gradient-free B 
modes [2, 3]. B modes are not generated at linear or- 
der in perturbation theory by the scalar perturbations 
which are the dominant source of CMB temperature and 
E-mode anisotropies. Because of this, B modes are of 
great interest as a clean probe of two more subtle sig- 
nals: (1) primordial tensor perturbations in the early 
Universe [4, 5], the measurement of which would provide 
a unique probe of the energy scale of inflation; and (2) 
gravitational lensing, which generates a distinctive non- 
Gaussian B-mode signal [6] that can be used to measure 
the projected mass distribution and constrain cosmolog- 
ical parameters such as the sum of neutrino masses (for 
a review, see [7]). 

Previous experiments have placed upper limits on the 
B-mode polarization anisotropy [8-11]. In this Letter we 
present the first detection of B modes sourced by grav- 
itational lensing, using first-season data from SPTpol, 
the polarization-sensitive receiver on the South Pole Tele- 
Scope. 

Gravitational lensing remaps the observed position of 
CMB anisotropies as ñ > fi4- Vo(ft), where ¢ is the CMB 
lensing potential [12]. This remapping mixes some of the 
(relatively) large E-mode signal into B. The induced B 
mode at Fourier wavevector lg is given to first order in 
$ as [13] 


pq) = [ ei | PGW G is iBT), 

(1) 
where the weight function W® specifies the mixing. In 
this Letter, we use measurements of E and 6$ to synthe- 
size an estimate for the lensing contribution, which we 
cross-correlate with measured B modes. Using maps of 
the cosmic infrared background measured by the SPIRE 
instrument onboard the Herschel space observatory to 
estimate ¢, and measurements of the E- and B-mode 
polarization from SP'Tpol, we detect the lensing signal 
at T.7c significance. 

CMB Data: The South Pole Telescope (SPT) [14] 
is a 10-meter telescope located at the geographic South 
Pole. Here we use data from SPTpol, a polarization- 
sensitive receiver installed on the telescope in January 
2012. SP'Tpol consists of two arrays of polarization- 
sensitive bolometers (PSBs): 1176 PSBs that observe at 
150 GHz, and 360 PSBs that observe at 95 GHz. The in- 
strument and its performance are described in Refs. [15- 
18]. The observation strategy, calibration, and data re- 
duction for SP Tpol data are similar to those used for the 
SPT-SZ survey, described in Ref. [19]. Here we briefly 
summarize the important points. 

We calibrate the PSB polarization sensitivities with 
observations of a ground-based thermal source behind a 
polarizing grid. This allows us to measure the polariza- 


tion angle of individual PSBs with « 2? statistical un- 
certainty and the average angle of all PSBs with « 0.1? 
statistical uncertainty. We estimate systematic uncer- 
tainty on the average angle to be « 1? (1.5?) at 150 GHz 
(95 GHz). 

Between March and November 2012, we used SPTpol 
to observe a 100 deg? region of low-foreground sky, be- 
tween 23h and 24h in right ascension and —50 and 
—60 deg in declination. We process the SP'Tpol data 
by “observations”, which are half-hour periods in which 
the telescope scans half of the field. Each observation 
is recorded as time-ordered data (TOD) from each PSB, 
in azimuthal scans separated by steps in elevation. For 
each scan, we apply a low-pass anti-aliasing filter as well 
as a high-pass 4th-order polynomial subtraction to re- 
move large-scale atmospheric fluctuations. This sup- 
presses modes along the scan direction, which we account 
for with a two-dimensional transfer function measured 
from simulations of the filtering process. 

In each observation, we drop PSBs with cuts based on 
noise level during the observation, response to elevation- 
dependent atmospheric power, and response to an in- 
ternal thermal calibration source. Typical observations 
include ~ 800 PSBs (~ 230 PSBs) at 150 GHz (95 GHz). 
We cut scans for PSBs with glitches (caused, for exam- 
ple, by cosmic ray hits). In typical 150 GHz (95 GHz) 
observations, we lose ~ 1% (~ 4%) of the data due to 
glitch removal. 

Data from each PSB are accumulated into maps of the 
I, Q, and U Stokes parameters using measured polar- 
ization angles and polarization efficiencies. We weight 
the TOD for each PSB in a scan by the inverse of 
the variance along the scan direction between 1 Hz and 
3Hz (1300 <1, < 3900 for the telescope scan speed of 
0.28 deg/s). We make maps using the oblique Lambert 
azimuthal equal-area projection [20] with square 2' x 2' 
pixels. This projection preserves area on the sky but in- 
troduces small distortions in angle; we account for these 
distortions by rotating the Q and U components to main- 
tain a consistent angular coordinate system across the 
map. For each observation we form a noise map from 
the difference of left- and right-going scans, cutting ob- 
servations which are outliers in metrics such as overall 
variance. This cut removes ^ 896 (7996) of the 150 GHz 
(95 GHz) data. Finally, we add the individual obser- 
vations together to produce full-season maps, with po- 
larization noise levels of approximately 10 uK-arcmin at 
150 GHz and 25 wK-arcmin at 95 GHz. 

Inaccuracy in PSB gain measurements can cause direct 
leakage of the CMB temperature into polarization, which 
we fit for using the cross-spectra of I with Q and U. We 
find « 296 leakage at both 150 GHz and 95 GHz, which 
we correct for by subtracting appropriate fractions of I 
from Q and U. We show below that this correction is 
unimportant for our final results. 

We calibrate the overall amplitude of the SPTpol maps 
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Left: Wiener-filtered E-mode polarization measured by SPTpol at 150GHz. Center: Wiener-filtered CMB lensing 


potential inferred from CIB fluctuations measured by Herschel at 500 um. Right: gravitational lensing B-mode estimate 
synthesized using Eq. (1). The lower left corner of each panel indicates the blue(-)/red(+) color scale. 


to better than 1% in temperature by cross correlating 
with SPT-SZ temperature maps over the same region 
of sky. The SPT-SZ maps are calibrated by comparing 
to the Planck surveyor 143 GHz maps [21] over the full 
2500 deg? SPT-SZ survey region. 

CIB Data: We use maps of the cosmic infrared back- 
ground (CIB) [22] obtained from the SPIRE instrument 
[23] onboard the Herschel space observatory [24] as a 
tracer of the CMB lensing potential ¢. The CIB has 
been established as a well-matched tracer of the lens- 
ing potential [22, 25, 26] and currently provides a higher 
signal-to-noise estimate of o than is available with CMB 
lens reconstruction. Its use in cross-correlation with the 
SPTpol data also makes our measurement less sensitive 
to instrumental systematic effects [27]. We focus on the 
Herschel 500 um map, which has the best overlap with 
the CMB lensing kernel [22]. 


Post-Map Analysis: We obtain Fourier-domain 
CMB temperature and polarization modes using a 
Wiener filter (for example, Ref. [28] and references 
therein), derived by maximizing the likelihood of the 
observed I, Q, and U maps as a function of the fields 
T(i), E(t), and B(I). The filter simultaneously decon- 
volves the two-dimensional transfer function due to the 
beam, TOD filtering, and map pixelization while down- 
weighting modes that are “noisy” due to either atmo- 
spheric fluctuations, extragalactic foreground power, or 
instrumental noise. We place a prior on the CMB auto- 
spectra, using the best-fit cosmological model given by 
Ref. [29]. We use a simple model for the extragalac- 
tic foreground power in temperature [19]. We use jack- 
knife difference maps to determine a combined atmo- 
sphere--instrument noise model, following Ref. [30]. We 
set the noise level to infinity for any pixels within 5' of 
sources detected at > 5e in [31]. We extend this mask to 
10' for all sources with flux greater than 50 mJy, as well as 
galaxy clusters detected using the Sunyaev-Zel’dovich ef- 


fect in Ref. [32]. These cuts remove approximately 5 deg? 
of the total 100deg? survey area. We remove spatial 
modes close to the scan direction with an £, < 400 cut, 
as well as all modes with | > 3000. For these cuts, our es- 
timated beam and filter map transfer functions are within 
2096 of unity for every unmasked mode (and accounted 
for in our analysis in any case). 

The Wiener filter naturally separates E and B con- 
tributions, although in principle this separation depends 
on the priors placed on their power spectra. To check 
that we have successfully separated E and B, we also 
form a simpler estimate using the yg formalism advo- 
cated in [33]. This uses numerical derivatives to estimate 
a field xp(x) which is proportional to B in harmonic 
space. This approach cleanly separates E and B, al- 
though it can be somewhat noisier due to mode-mixing 
induced by point source masking. We therefore do not 
mask point sources when applying the yg estimator. 

We obtain Wiener-filtered estimates QCIP of the lens- 
ing potential from the Herschel 500 jum maps by apply- 
ing an apodized mask, Fourier transforming, and then 
multiplying by QUT Cee er). We limit our 
analysis to modes l > 150 of the CIB maps. We model 
the power spectrum of the CIB following Ref. [34], with 
Cf1P-CIB. — 3500(1/3000)-L25Jy?/sr. We model the 
cross-spectrum Cf |"? between the CIB fluctuations and 
the lensing potential using the SSED model of Ref. [35], 
which places the peak of the CIB emissivity at redshift 
Ze = 2 with a broad redshift kernel of width o, = 2. 
We choose a linear bias parameter for this model to 
agree with the results of Refs. [22, 26]. More realistic 
multi-frequency CIB models are available (for example, 
Ref. [36]); however, we only require a reasonable tem- 
plate. The detection significance is independent of errors 
in the amplitude of the assumed CIB-¢ correlation. 

Results: In Fig. 1, we plot Wiener-filtered estimates 
E99 and ¢°'® using the CMB measured by SPT pol at 
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FIG. 2: Black, center bars: cross correlation of the lensing 


B modes measured by SPTpol at 150 GHz with lensing B 
modes inferred from CIB fluctuations measured by Herschel 
and E modes measured by SPTpol at 150 GHz; as shown in 
Fig. 1. Green, left-offset bars: same as black, but using E 
modes measured at 95 GHz, testing both foreground contam- 
ination and instrumental systematics. Orange, right-offset 
bars: same as black, but with B modes obtained using the 
xB procedure described in the text rather than our fiducial 
Wiener filter. Gray bars: curl-mode null test as described in 
the text. Dashed black curve: lensing B-mode power spec- 
trum in the fiducial cosmological model. 


150 GHz and the CIB fluctuations traced by Herschel. 
In addition, we plot our estimate of the lensing B modes 
obtained by applying Eq. (1) to these measurements. In 
Fig. 2 we show the cross-spectrum between this lensing 
B-mode estimate and the B modes measured directly 
by SPTpol. The data points are a good fit to the ex- 
pected cross-correlation, with a x?/dof of 3.5/4 and a 
corresponding probability-to-exceed (PTE) of 48%. We 
determine the uncertainty and normalization of the cross- 
spectrum estimate using an ensemble of simulated, lensed 
CMB-+noise maps and simulated Herschel maps. We 
obtain comparable uncertainties if we replace any of the 
three fields involved in this procedure with observed data 
rather than a simulation, and the normalization we de- 
termine for each bin is within 15% of an analytical pre- 
diction based on approximating the Wiener filtering pro- 
cedure as diagonal in Fourier space. 


In addition to the cross-correlation E$ x B, it is also 
interesting to take a “lensing perspective” and rear- 
range the fields to measure the correlation EB x $. In 
this approach, we perform a quadratic “EB” lens re- 
construction [13] to estimate the lensing potential PEB, 
which we then cross-correlate with CIB fluctuations. The 
observed cross-spectrum can be compared to previous 
temperature-based lens reconstruction results [22, 26]. 
This cross-correlation is plotted in Fig. 3. Again, the 
shape of the cross-correlation which we observe is in good 
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FIG. 3: “Lensing view” of the EBọ correlation plotted in 
Fig. 2, in which we cross-correlate an EB lens reconstruc- 
tion from SPTpol data with CIB intensity fluctuations mea- 
sured by Herschel. Left green, center black, and right or- 
ange bars are as described in Fig. 2. Previous analyses using 
temperature-based lens reconstruction from Planck [26] and 
SPT-SZ [22] are shown with boxes. The results of Ref. [26] are 
at a nominal wavelength of 550 um, which we scale to 500 jum 


with a factor of 1.22 [37]. The dashed black curve gives our 


fiducial model for a "B® as described in the text. 


agreement with the fiducial model, with a y?/dof of 2.2/4 
and a P'TE of 7096. 

Both the E$ x B and EB x ¢ cross-spectra discussed 
above are probing the three-point correlation function 
(or bispectrum) between E, B, and ¢ that is induced by 
lensing. We assess the overall significance of the mea- 
surement by constructing a minimum-variance estimator 
for the amplitude A of this bispectrum, normalized to 
have a value of unity for the fiducial cosmology+CIB 
model (analogous to the analyses of Refs. [38, 39] for 
the T'T'ó bispectrum). This estimator can be written as 
a weighted sum over either of the two cross-spectra al- 
ready discussed. Use of A removes an arbitrary choice 
between the lensing or B-mode perspectives, as both are 
simply collapsed faces of the EB bispectrum. Relative 
to our fiducial model, we measure a bispectrum ampli- 
tude A = 1.092 + 0.141, non-zero at approximately 7.70. 

We have tested that this result is insensitive to analy- 
sis choices. Replacement of the B modes obtained using 
the baseline Wiener filter with those determined using 
the xp estimator causes a shift of 0.20. Our standard 
B-mode estimate incorporates a mask to exclude bright 
point sources, while the xp estimate does not. The good 
agreement between them indicates the insensitivity of po- 
larization lensing measurements to point-source contam- 
ination. If we change the scan direction cut from lẹ < 400 
to 200 or 600, the measured amplitude shifts are less than 
1.2c, consistent with the root-mean-squared (rms) shifts 
seen in simulations. If we repeat the analysis without 


correcting for I + Q, U leakage, the measured amplitude 
shifts by less than 0.1c. A similar shift is found if we ro- 
tate the map polarization vectors by one degree to mimic 
the E — B leakage which would be induced by an error 
in the average PSB angle. 

We have produced estimates of Bleus using alterna- 
tive estimators of E. When we replace the E modes 
measured at 150 GHz with those measured at 95 GHz, 
we measure an amplitude A=1.225+ 0.164, indicating 
a lack of significant foreground contamination or spuri- 
ous correlations between 150 GHz E and B modes. We 
have also estimated the bispectrum amplitude using the 
CMB temperature as a tracer of E (exploiting the CIE 
correlation, as in Refs. [40, 41]). We obtain an ampli- 
tude of Â = 1.374 + 0.427, consistent with unity. The 
shifts between these values of Â and our fiducial one are 
consistent with the rms scatter seen in simulations. 

Our amplitude estimate also passes several null 
tests. Using E modes obtained from the difference 
map of left-going and right-going scans, we measure 
A = 0.049 + 0.037. Replacing instead the B modes with 
a difference map, we obtain A = —0.028 + 0.119. We also 
estimate the amplitude of the B-mode power sourced by 
curl-type lensing modes. These generate deflections as 
the curl of a scalar potential rather than a gradient [42]. 
Curl-type lensing modes are expected to be negligibly 
small in the fiducial cosmological model [43], but could 
be sourced by instrumental effects. If we take the CIB 
as a tracer of curl-type lensing modes and measure the 
lensing B-mode power in cross-correlation analogously, 
we obtain the gray errorbars in Fig. 2, which are consis- 
tent with the expected value of zero, having a x?/dof of 
2.9/4 (PTE = 57%). 

As a final consistency test, we also replace the CIB- 
based estimate of @ with quadratic lens reconstruc- 
tions based on combinations TT, TE, EE, and EB 
of the CMB fluctuations [13, 27, 43, 44]. These lens- 
ing estimators probe the CMB lensing potential di- 
rectly, and are not limited by the imperfect overlap 
of the CIB redshift kernel with that of lensing. The 
temperature-only lensing estimator has already been 
used to make precise measurements of the lensing po- 
tential power spectrum [30, 45-47]. We now have the 
ability to extend this analysis to polarization. We have 
not yet performed a thorough characterization of these 
estimators or their sensitivity to analysis choices and 
possible systematic effects in the SPTpol data; how- 
ever, comparing the measured amplitudes with our fidu- 
cial result using statistical error bars is still a use- 
ful consistency test. We obtain amplitude estimates 
in reasonable agreement with the fiducial cosmological 
model, with A?’ =0.675 + 0.194, AT” =0.921 + 0.294, 
AE —0.860 + 0.267, and AE? = 0.960 + 0.386. These 
estimators probe the four-point correlations (or trispec- 
tra) of the lensed CMB, which can be non-zero even 
for Gaussian fluctuations. We have subtracted a “bias- 


hardened” estimate of the Gaussian contribution [48] for 
each estimator to correct for this. 


Conclusions: B-mode polarization is a promising av- 
enue for measuring both the CMB lensing potential ¢ 
and primordial gravitational waves. In this Letter, we 
have presented the first detection of B modes produced 
by gravitational lensing. Beyond this detection, the lens- 
ing B-mode map which we have synthesized can be sub- 
tracted from the observed B modes in a process of “de- 
lensing”, reducing the effective noise level for the mea- 
surement of primordial B modes [7, 49-51]. The work 
presented here is a first step in the eventual exploitation 
of CMB B-mode polarization as a probe of both structure 
formation and the inflationary epoch. 


Acknowledgements: The SPT is supported by 
the National Science Foundation through grant ANT- 
0638937, with partial support provided by NSF grant 
PHY-1125897. Support for the development and con- 
struction of SPTpol were provided by the Gordon and 
Betty Moore Foundation through Grant GBMF 947 to 
the University of Chicago, a gift from the Kavli Foun- 
dation, and NSF award 0959620. Herschel is an ESA 
space observatory with science instruments provided by 
European-led Principal Investigator consortia and with 
important participation from NASA. This research used 
resources of the National Energy Research Scientific 
Computing Center, which is supported by the Office of 
Science of the U.S. Department of Energy under Con- 
tract No. DE-AC02-05CH11231. It also used resources 
of the CLUMEQ supercomputing consortium, part of the 
Compute Canada network. Research at Argonne Na- 
tional Laboratory and use of the Center for Nanoscale 
Materials are supported by the Office of Science of the 
U.S. Department of Energy under Contract DE-AC02- 
06CH11357. The McGill group acknowledges funding 
from the National Sciences and Engineering Research 
Council of Canada, Canada Research Chairs program, 
and the Canadian Institute for Advanced Research. The 
C.U. Boulder group acknowledges support from NSF 
AST-0956135. We thank P. Hargrave at Cardiff Univer- 
sity for anti-reflection coating the SPT pol lens, A. Dates- 
man for his work on TES detectors at Argonne, R. Divan 
for microfabrication support at Argonne, and the mem- 
bers of the Truce collaboration for their efforts in the de- 
sign of the 150 GHz polarization detectors, in particular 
D. Becker, J. Britton, M.D. Niemack, and K.W. Yoon at 
NIST. We thank M. Lueker, T. Plagge, Z. Staniszewski, 
E. Shirokoff, H. Spieler and R. Williamson for their con- 
siderable contributions to the SPT program. DH was 
supported by the Lorne Trottier Chair in Astrophysics 
and Cosmology at McGill as well as a CITA National 
Fellowship. RK acknowledges support from NASA Hub- 
ble Fellowship Grant HF-51275. 


[16] 


[17] 


[18] 


[19] 


[20] 


[21] 
[22] 


[23] 


W. Hu and S. Dodelson, Ann.Rev.Astron.Astrophys. 40, 
171 (2002), astro-ph/0110414. 

M. Kamionkowski, A. Kosowsky, and A. Stebbins, 
Phys.Rev. D55, 7368 (1997), astro-ph/9611125. 

M. Zaldarriaga and U. Seljak, Phys.Rev. D55, 1830 
(1997), astro-ph/9609170. 

M. Kamionkowski, A. Kosowsky, and A. Stebbins, 
Phys.Rev.Lett. 78, 2058 (1997), astro-ph/9609132. 

U. Seljak and M. Zaldarriaga, Phys.Rev.Lett. 78, 2054 
(1997), astro-ph/9609169. 

M. Zaldarriaga and U. Seljak, Phys.Rev. D58, 023003 
(1998), astro-ph/9803150. 

K. M. Smith, A. Cooray, S. Das, O. Dore, D. Hanson, 
et al., AIP Conf.Proc. 1141, 121 (2009), 0811.3916. 

H. Chiang, P. Ade, D. Barkats, J. Battle, E. Bierman, 
et al., Astrophys.J. 711, 1123 (2010), 0906.1181. 

M. Brown et al. (QUaD collaboration), Astrophys.J. 705, 
978 (2009), 0906.1003. 

C. Bennett et al. (WMAP Collaboration) 
1212.5225. 

D. Araujo et al. (QUIET Collaboration), Astrophys.J. 
760, 145 (2012), 1207.5034. 

A. Lewis and A. Challinor, Phys.Rept. 429, 1 (2006), 
astro-ph/0601594. 

W. Hu and T. Okamoto, Astrophys.J. 574, 566 (2002), 
astro-ph/0111606. 

J. Carlstrom, P. Ade, K. Aird, B. Benson, L. Bleem, 
et al., Publ.Astron.Soc.Pac. 123, 568 (2011), 0907.4445. 
E. M. George, P. Ade, K. A. Aird, J. E. Auster- 
mann, J. A. Beall, D. Becker, A. Bender, B. A. Ben- 
son, L. E. Bleem, J. Britton, et al., in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference 
Series (2012), vol. 8452, 1210.4971. 

J. E. Austermann, K. A. Aird, J. A. Beall, D. Becker, 
A. Bender, B. A. Benson, L. E. Bleem, J. Britton, J. E. 
Carlstrom, C. L. Chang, et al., in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference 
Series (2012), vol. 8452, 1210.4970. 

J. T. Sayre, P. Ade, K. A. Aird, J. E. Austermann, J. A. 
Beall, D. Becker, B. A. Benson, L. E. Bleem, J. Britton, 
J. E. Carlstrom, et al., in Society of Photo- Optical Instru- 
mentation Engineers (SPIE) Conference Series (2012), 
vol. 8452, 1210.4968. 

J. W. Henning, P. Ade, K. A. Aird, J. E. Austermann, 
J. A. Beall, D. Becker, B. A. Benson, L. E. Bleem, 
J. Britton, J. E. Carlstrom, et al., in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference 
Series (2012), vol. 8452, 1210.4969. 

K. T. Story, C. L. Reichardt, Z. Hou, R. Keisler, K. A. 
Aird, B. A. Benson, L. E. Bleem, J. E. Carlstrom, 
C. L. Chang, H. Cho, et al., Submitted to ApJ (2012), 
1210:7231. 

J. P. Snyder, Map projections-a working manual (U.S. 
Geological Survey, Washington, 1987). 

P. Ade et al. (Planck Collaboration) (2013), 1303.5062. 
G. Holder, M. Viero, O. Zahn, K. Aird, B. Benson, et al., 
Astrophys.J. 771, L16 (2013), 1303.5048. 

M. Griffin, A. Abergel, A. Abreu, P. Ade, P. André, et al., 
Astron.Astrophys. 518, L3 (2010), 1005.5123. 


(2012), 


24 


25 


26 
27 


28 


29 
30 


31 


32 


G. Pilbratt, J. Riedinger, T. Passvogel, G. Crone, 
D. Doyle, et al., Astron.Astrophys. 518, L1 (2010), 
1005.5331. 

Y.-S. Song, A. Cooray, L. Knox, and M. Zaldarriaga, 
Astrophys.J. 590, 664 (2003), astro-ph/0209001. 

P. Ade et al. (Planck Collaboration) (2013), 1303.5078. 
K. Benabed, F. Bernardeau, and L. Van Waerbeke, 
Phys.Rev. D63, 043501 (2001), astro-ph/0003038. 

F. Elsner and B. D. Wandelt, Astron.Astrophys. 549, 
A111 (2013), 1210.4931. 

P. Ade et al. (Planck Collaboration) (2013), 1303.5076. 
A. van Engelen, R. Keisler, O. Zahn, K. Aird, B. Benson, 
et al., Astrophys.J. 756, 142 (2012), 1202.0546. 

L. M. Mocanu, T. M. Crawford, J. D. Vieira, K. A. 
Aird, M. Aravena, J. E. Austermann, B. A. Benson, 
M. Béthermin, L. E. Bleem, M. Bothwell, et al., ArXiv 
e-prints (2013), 1306.3470. 

K. Vanderlinde, T. M. Crawford, T. de Haan, J. P. Dud- 
ley, L. Shaw, P. A. R. Ade, K. A. Aird, B. A. Benson, 
L. E. Bleem, M. Brodwin, et al., ApJ 722, 1180 (2010), 
1003.0003. 

K. M. Smith and M. Zaldarriaga, Phys.Rev. D76, 043001 
2007), astro-ph/0610059. 

G. E. Addison, J. Dunkley, A. Hajian, M. Viero, J. R. 
Bond, et al., Astrophys.J. 752, 120 (2012), 1108.4614. 
N. R. Hall, L. Knox, C. L. Reichardt, P. A. R. Ade, K. A. 
Aird, B. A. Benson, L. E. Bleem, J. E. Carlstrom, C. L. 
Chang, H. Cho, et al., ApJ 718, 632 (2010), 0912.4315. 
M. Béthermin, L. Wang, O. Doré, G. Lagache, M. Sar- 
gent, et al. (2013), 1304.3936. 

R. Gispert, G. Lagache, and J. Puget, A&A 360, 1 
2000), astro-ph/0005554. 

K. M. Smith, O. Zahn, and O. Dore, Phys.Rev. D76, 
043510 (2007), 0705.3980. 

C. M. Hirata, S. Ho, N. Padmanabhan, U. Seljak, 
and N. A. Bahcall, Phys.Rev. D78, 043520 (2008), 
0801.0644. 

A. H. Jaffe, New Astron.Rev. 47, 1001 (2003), 
ph/0307118. 

O. Dore, G. P. Holder, and A. Loeb, Astrophys.J. 612, 
81 (2004), astro-ph/0309281. 

A. Cooray, M. Kamionkowski, and R. R. Caldwell, 
Phys.Rev. D71, 123527 (2005), astro-ph/0503002. 

C. M. Hirata and U. Seljak, Phys.Rev. D68, 083002 
(2003), astro-ph/0306354. 

J. Guzik, U. Seljak, and M. Zaldarriaga, Phys.Rev. D62, 
043517 (2000), astro-ph/9912505. 

S. Das, B. D. Sherwin, P. Aguirre, J. W. Appel, 
J. R. Bond, et al., Phys.Rev.Lett. 107, 021301 (2011), 
1103.2124. 

S. Das, T. Louis, M. R. Nolta, G. E. Addison, E. S. 
Battistelli, et al. (2013), 1301.1037. 

P. Ade et al. (Planck Collaboration) (2013), 1303.5077. 
T. Namikawa, D. Hanson, and R. Takahashi, MNRAS 
431, 609 (2013), 1209.0091. 

L. Knox and Y.-S. Song, Phys.Rev.Lett. 89, 011303 
(2002), astro-ph/0202286. 

M. Kesden, A. Cooray, and M. Kamionkowski, 
Phys.Rev.Lett. 89, 011304 (2002), astro-ph/0202434. 
U. Seljak and C. M. Hirata, Phys.Rev. D69, 043005 
(2004), astro-ph/0310163. 


astro- 


